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This report deals with scattering of electromagnetic waves from an 
infinite cylinder of uniform but complex dielectric constant. Its main 
purpose is to make available a set of curves showing the dependence of 
transmission and reflection coefficients and their phases on plasma fre- 
quency and cylinder radius. Two different collision frequencies are 
used. All quantities have been normalized to eliminate the dependence 
on the operating frequency and the distance between the transmitting and 
receiving horns. The incident wave is assumed to be a plane wave travel- 
ing in a direction normal to the cylinder axis. Both the transverse mag- 
netic. case and the transverse electric case are: considered. 
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The solution of a scattering problem involves three steps 
1) Determination of the electric field E from the equation 


Peer? e Be UV-B) (2) 
or of the magnetic field H from the equation 


Fut e H= (VE)XS- (2) 


where k, =w CH, is the free space propagation constant and ¢' is the 
camplex dielectric constant*. Which of the two equations 1 or 2 is 
_#olved depends on which is simpler for a given geometry and incident 
rediation. Only one equation has to be solved because once E has been 
found H can be determined with the help of Maxwell's equations end simi- 
larly once H has been found E can be determined. 

2) Satisfaction of the radiation condition at iafinity. 

3) Satisfaction of the boundary conditions at the interface between 
regions consisting of different media. 


* See appendixi 
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Coordinate System Cross Section of p a cylinder 
with axis coinciding with the z-axis 


Outside the cylinder «' = 1 and inside the cylinder ¢’ is taken to be a 
constant given by 
2 2 (1.45) 
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where 2, = -B is the ratio of the plasma frequency f, to f the frequency 
of the incident wave, 
only = 345 is the ratio of y the electron collision frequency to xf 
the circular frequency of the incident wave. 

For the transverse magnetic case the incident electric field is 
parallel to the z-axis. In this case equation 1 is simpler than equation 
2 and reduces to the scalar equation 
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Neglecting the time depedence et the incident field is given by 


eo 
n - 
Eo x i J, (xr) cosne 
n=o 


where i, is a unit vector in the s direction 
r sz = 
ay 
Ls the free space wave length of the incident radiation 
1 nz0o 

es 2 when | >1 

Let the scattered field outside the cylinder be represented by 
o 


E = 18, ). re a(2) (27) cosno 
m0 


and the total field inside the cylinder by 


o 
= 18, y Cet dF (ant /et ) cosne 
n=o 

The choice of the Bessel function in BY insures that the solution will 
be finite at the origin and that of the Hankel function of the first kind 
in £° insures that the scattered field at large distances from the cylin- 
der axis will behave as an outgoing cylindrical wave and the radiation 
condition will be satisfied. 

The constants A, and C. can be determined from the boundary condi- 


tions. To obtain reflection and transmission coefficients and their 


wey 


Phases the scattered field outside the cylinder must be known and hence 
A, must be found. Using Mexwe11's equation we find 


1 iE n ayn ear) 
Hy i » = 3(aaF) cos nO 
o (1) - 
i = = ». Aen ng — cosnd 
lo 
, Ie J (anrer) 
Hy = 7 sy aw ar haa 


where 1% “fee is the impedance of free space. 
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The boundary conditions at rae require that 
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Substituting the series expressions for the field components and letting 


.-« 2 we obtain 
Ay 


AB.) (ont) + J, (2nd) = co. (and 4K") 


an?) (ant) ag. (oxi) a5,(axi er) 
Bent)” tens) Sea ae) 


From these using the expression for the derivative of a Bessel function 
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where 


&" Jy (2x8 A”) N(2x8) - J (2x0 ve” ) W4 (2xa) 
tn” er ty (2nd 4) 5 (2nd) - J (2nd 7) 7, (Ona) 
Similarly for the transverse electric case the incident magnetic field 


is parallel to the z-axis, and neglecting time dependence is given by 


E e Ee 
i o .12xrcos0 fe} n - 
aa eo a oa y, ei” J (2xt) cosnd 


nso 


proceeding as before let the scattered magnetic field outside the cylin- 
der be given by 


E eo 

5. n .(1) - 

E, os Bet H, (2ar) cosnd 
n2zo 


and the total magnetic field inside the cylinder by 


E 
t ‘oO 7 n é 
a i, fe De ,i Jy(2xr 7") cosnd 


using Maxwell’s equations we find 


n 2, (2x7) 


1 
E wz «AE 3 ei ———— + ¢cosnd 
e ° n - 
feo 0 (2ar) 


g° « “iE, Bet on) (ant) cosne 
n=o 0 (2ar) 


= OJ (2xrj 
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For the transverse electric case the boundary conditions which mst be 


satisfied at r=a give 


a5,(28) ag Mare) a,x Ye”) 
+ —_——- «= ———— 
3(2ma) nO (28) nd (20a Je”) 


J(2n8) + BH? )(ana) = 4" DJ (288 67) 


From which 


where 


Jy (2nd YEN (2n8)~ 2S (2nd Ve IN, (2n8)4B ef el 
b - : : seat Te 
n Jy (2nd Ye )s (2ne)~ 7S, (2na fe” 7, (2x8) Be 


J (2x8 Ye” )N (2x8) 
n(2xa /e”)J, (2x8) 


The program makes use of an JEM library subroutine to calculate the 
Bessel and Newnann functions for real arguments and of the series expres- 
sion together with the recursion forma for the Bessel functions to 
determine the Bessel functions for complex arguments. Care has been 
taken to avoid overflow and round off troubles in the calculation of A, 
and Be 

Let 
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EB at a large distance from the axis of the cylinder 


eA the scattered electric field for the transverse magnetic case /0 = 0 
Onux 


the scattered magnetic field for the transverse electric case = O 


GS 
gE at a large distance from the axis of the cylinder =x 


-_ 
E'). distance from the axis of the cylinder =x 


the incident magnetic field at a large = 0 


z distance from the axis of the cylinder =x 
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the incident electric field at a large 020 
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We define the following eight quantities 


2 e J 
1 = 2x°F 5 = : €n(An tno) ; + ». €n(AjoAn) 
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PIM = phase “(2 = tan B= 
: D, fala + Ay) 
=f Bo Bn 


- ) eg(-2 a2 
PRM = phase of (3 +x/h- har = ten um 
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these are the transmission and reflection coefficients and phases 
appearing in the graphs. 


I. Dielectric Constant 


Lut 


Assuming a time dependence e Maxwell's equations are 


VXE = iu 
V¥H = -lue E+ J 


where €. and H, are the dielectric constant and permeability of free 


space respectively. For a single electron 


v being the electron collision frequency and m and e the mass and charge 
of the electron. 
Hence 


ieE (1-i¥) 
x = 


we 
<s 
« 

el< 


mu(l + 9°) 


If J = Ne x, where N is the electron density, then 


YX H = -iwe E (1 - 
= = me) (1+y 


or 


VXH = -iwe e' E 


where 


2 a 
Pu er its Ne j-tv) 


mew (1+¥") 
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This is the effective complex dielectric constant. 


The quantity t> = ee is usually referred to as the plasma 
em 
{e) . 


frequency of the medium. 


rf 
If #, = x, where f is the frequency of the incident radiation 
then 
2? (1-45) 


This is the form in which the dielectric constant is used in the progran. 


i and v are input values and e’ is calculated. 
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Figures 1 - 2h 
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Figure 1 
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Figure 2 
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Figure 4 
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Figure 5 
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Figure 7 
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Figure 9 
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APPENDIX IT 


The PRM and PRE curves possess a few interesting characteristics 
which will be briefly discussed here. It appears that for very low elec- 
tron densities, i.e. for very small ae the phase of the reflection coef- 
ficient for both the transverse magnetic case and the transverse electric 
case becomes independent of the electron density. Given a normalized 
collision frequency vy and any normalized radius a PRM always takes one 
of two values which differ by 2, and the same is true for PRE. 

Why this is so can easily be seen from the equations arising from 
the solution of the boundary value problem. Assuming a to be small 


we may write 


# 2 (1.43) 


ei = 1] - = 1-5 
145° 
and 
7 8 
Ve" = 1-5 
where 
2 = (1-43) 
8 = 
iets 


Discarding all terms involving powers of 8 higher then the 2nd we obtain 


re [BAS K (2nd) 
- tan v + tan aa Sear eeuenena 
~J, (4a) 


where 
co 2 
iat jks n Sie pean: cele 
K(2na) = Ona ) (-1) J, (ana) 7g (28) ~ J, a) 


( PRM + PRE but PRM = PRE + x because we have defined the reflection coef-~ 
ficient in terms of the scattered electric field in the transverse mag- 
netic case and in terms of the scattered magnetic field in the trans- 


verse electric case +) 


This expression shows at once that for small ae 


1 


PRM = PRE + = - = - tan” ¥ when a is such that J, (4xa) > 0 


njA 


: 
: 


~ tan”! > when & is such thet J, (4na) << 0 


PRM = PRE +x = - tan”! } when & is such that J, (42a) = O and K(2na)< 0 
1 


-tan> 


2 
z 


v when a is such that Jy (4xa) = O and K(2nxa)> 0 


also 


2 


Hence for small a: the slope of a curve representing PRM or PRE as a 
function of a depends for a given & on both J, (4a) and K(2xa). Suppose 
a is varied by small increments until J, (4a) changes sign. Since K(2na) 
will not change sign at the same time as J, (nt. ), PRM or PRE will change 
by 2, and the slope of a curve representing PRM or PRE as a function of 
ee will change sign. If a is varied further K(2xa) will eventually 


change sign too and the slope will change sign again. For a such that 


A-2 


J, (4xa) = 0 the slope of a curve representing PRM or PRE as a function 
of a will be + » or -o depending on whether K(2na) is positive or 
negative. Figure 78 shows in some detail the behavior of PRM for small 


a) a x - 
fo when 4 is near a value a, where J, (4na, ) = 0, 
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